Nucleon form factors and a nonpointlike diquark 
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Nucleon form factors are calculated on q 2 £ [0, 3] GeV 2 
using an Ansatz for the nucleon's Fadde'ev amplitude moti- 
vated by quark-diquark solutions of the relativistic Fadde'ev 
equation. Only the scalar diquark is retained, and it and the 
quark are confined. A good description of the data requires 
a nonpointlike diquark correlation with an electromagnetic 
radius of 0.8 r n . The composite, nonpointlike nature of the 
diquark is crucial. It provides for diquark-breakup terms that 
are of greater importance than the diquark photon absorption 
contribution. 

Mesons present a two-body problem, and the Dyson- 
Schwinger equations (DSEs) have been widely used in 
the calculation of their properties and interactions GL^] . 
Many studies have focused on electromagnetic processes; 
such as the form factors of light pseudoscalar and 
vector mesons and the 7*7r° — > 7 |5|-[|, 7*7r — » p |t]] 
and 77T* — ► irir |B| transition form factors, all of which are 
accessible at TJNAF. These studies provide a foundation 
for the exploration of nucleons, which is fundamentally a 
three-body problem. 

The nucleon's bound state amplitude can be obtained 
from a relativistic Fadde'ev equation |[(|. Its analysis 
may be simplified by using the feature that ladder-like 
dressed-gluon exchange between quarks is attractive in 
the colour antitriplet channel. Then, in what is an ana- 
logue of the rainbow-ladder truncation for mesons, the 
Fadde'ev equation can be reduced to a sum of three cou- 
pled equations, in which the primary dynamical content 
is dressed-gluon exchange generating a correlation be- 
tween two quarks and the iterated exchange of roles be- 
tween the dormant and diquark-participant quarks. Fol- 
lowing this approach, the diquark correlation is repre- 
sented by the solution of an homogeneous Bethe-Salpeter 
equation in the dressed-ladder truncation and hence its 
contribution to the quark-quark scattering matrix, M. qq , 
is that of an asymptotic bound state; i.e., it contributes 
a simple pole. That is an artefact of the ladder trun- 
cationjO and complicates solving the Fadde'ev equa- 
tion |12| by introducing spurious free-particle singulari- 
ties in the kernel. 

Studies of DSE-models [QJ|] suggest that confinement 
can be realised via the absence of a Lehmann repre- 
sentation for coloured Green functions, and have led 
to a phenomenologically efficacious parametrisation of 
the dressed-quark Schwinger function ||. A similar 
parametrisation of the diquark contribution to M. qq , ad- 
vocated in Ref . (Of , has been used to good effect in solv- 



ing the Fadde'ev equation [jl4|. We use such representa- 
tions herein. 

The nucleon-photon current is[] 



J IM {P',P)=ieu(P')A fi (q,P)u(P). 



(1) 



where the spinors satisfy: 7 ■ Pu(P) = iMu(P), u(P) 7 • 
P = iMu(P), with M = 0.94 GeV the nucleon mass, and 
q = (P' — P). The complete specification of a fermion- 
vector-boson vertex requires twelve independent scalar 
functions: 

i\p{q, P) = /1 + iOp V q v f 2 + R a h + »7 • RRfi. fi 
+i(7 V p R u q v Rp h + il^lv Spvpcr q P Rcr h + ■ ■ • I ( 2 ) 

where /, - Mq 2 , q ■ P, P 2 ), R = (P' + P) and q ■ R = 
for elastic scattering. However, using the definition of 
the nucleon spinors, (pJ) can be written 



(3) 



ieu(P') (7,, W) + <V Q« W) J u(P) , 

where the Dirac and Pauli form factors are 

Ft = (4) 
fx + 2M/3-4 M 2 U-2M q 2 f B - q 2 f 6 , 
F2 = (5) 
2M/ 2 - 2M/ 3 + 4M 2 / 4 + 2 M h -AM 2 f 6 , 

in terms of which one has the electric and magnetic form 
factors: 



G E (q 2 )^F 1 (q 2 )-^_F 2 (q 2 ) 
G M (q 2 ) = F 1 (q 2 ) + F 2 (q 2 ). 



(6) 
(7) 



To calculate these form factors we represent the nu- 
cleon as a three-quark bound state involving a diquark 
correlation, and require the photon to probe the di- 
quark's internal structure. Antisymmetrisation ensures 
there is an exchange of roles between the dormant and 



1 In our Euclidean formulation: p-q = Pi<7i, {jp, 7^} 

■ = *[ln,'Yu\ > and tro [757^7,^7^] 



2 5p„, 7+ 



7m. 



-Atpvpa, £1234 = 1. 



1 



diquark-participant quarks and this gives rise to diquark 
"breakup" contributions. We describe the propagation 
of the dressed-quarks and diquark correlation by confin- 
ing parametrisations and hence pinch singularities asso- 
ciated with quark production thresholds are absent. Our 
calculation is kindred to many studies of meson proper- 
ties 

We write the Fadde'ev amplitude of the nucleon as |lq| 
^ T a (p 1 ,a 1 ,T 1 ;p 2 ,a 2 ,T 2 ;p 3 ,a 3 ,T 3 ) = (8) 

3 12 
£c±C2C3 & &CX.OC3 Ql Q>2 

where s ClC2C3 effects a singlet coupling of the quarks' 
colour indices, (pi, 014, t % ) denote the momentum and the 
Dirac and isospin indices for the i-th quark constituent, a 
and r are these indices for the nucleon itself, ip{(,\, ^2) is a 
Bethe-Salpeter-like amplitude characterising the relative- 
momentum dependence of the correlation between di- 
quark and quark, A (if) describes the propagation char- 
acteristics of the diquark, and 

r£«(pi>K0 = (Ci l5 ) aia2 (iT 2 y^ r( Pl , P2 ) (9) 

represents the momentum-dependence, and spin and 
isospin character of the diquark correlation; i.e., it corre- 
sponds to a diquark Bethe-Salpeter amplitude. 

With this form of we retain in M. qq only the contri- 
bution of the scalar diquark, which has the largest cor- 
relation length p3[ : A + := l/mo+ = 0.27 fm. For all 
(wd)-correlations with J p ^ 1 + , X u( j < 0.5 A +. The 
axial- vector correlation is different: Ai+ = 0.78 Ao+, and 
it is quantitatively important in the calculation of baryon 
masses (< 30%) [Q. Hence we anticipate that neglect- 
ing the 1 + correlation will prove the primary defect of 
our Ansatz. However, it is an helpful expedient in this 
exploratory calculation, which is made complicated by 
our desire to elucidate the effect of the diquarks' internal 
structure. 

Our impulse approximation to the nucleon form factor 
is depicted in Fig. |l|. Enumerating from top to bottom, 
the diagrams represent 

Al( q ,P) = 3j J 0z (10) 

$(K,pa + q)A{K)^(K,p 3 ) QfA«( P3 + g,pa) , 

witlfl K = V P+l,p 3 = {l- Tl)P-i, p 2 =K/2-k,Q F = 
diag(2/3, -1/3), A«(fci, fc 2 ) = S(k 1 )T^k 1 ,k 2 ) S(k 2 ), 

Al(q 7 P)=6j 1 0z 1 s n(p 1 +q,p 2 ,p 3 ) (11) 
xn(p 1 ,p 2 ,p 3 ) tr D [A^(pi + q,Pi)S(p 2 )] S(p 3 ) \I F , 



2 rj describes the partitioning of the nucleoli's total momen- 
tum: P = pi + p2 + P3 , between the diquark and quark, a 
necessary feature of a covariant treatment. 



which contributes equally to the proton and neutron and 
contains the diquark electromagnetic form factor, with 

6 — ^ciC2C3^ciC2C3 and 

tt(Pl,P2,P3) = i>(Pi +P2,Ps) A(pi +p 2 )T(p 1 ,p 2 ) , (12) 

kl(q,P)=& J^^n{p l+ q,p 3 ,p 2 ) (13) 
xn( Pl , P2 , P3 ) S{ P2 ) (iT 2 ) T Q F {iT 2 ) A«(pi,pi + q) S(p 3 ) , 

Afa,P) = 6j$jt$tn(pi,Ps,P2 + q) (W) 
xn(p 1 ,p 2 ,p 3 ) Q F A% (P2 + q,P2) S(pi) S{p 3 ) , 
A«(«,P) = 6 J J^ T 0zn( Pl ,p 3 + q ,p 2 ) (15) 
xn(p 1 ,p 2 ,p 3 ) S(p 2 ) 5(pi) <3fA m (p 3 + q,p 3 ) . 




FIG. 1. Our impulse approximation to the electroma gne tic 
current requires the calculation of five contributions, (not) - 
(p"5|). ip(li,l2) in (^); F: Bethe-Salpeter-like diquark am- 
plitude in (^); solid line: S(q), quark propagator in (|l7|); 
dotted line: A (if), diquark propagator in (^). The lowest 
three diagrams, which describe the interchange between the 
dormant quark and the diquark participants, effect the an- 
tisymmetrisation of the nucleon's Fadde'ev amplitude. Cur- 
rent conservation follows because the photon-quark vertex is 
dressed, given in (jj^). 

The nucleon-photon vertex is 

5 

^(q,P)=Al(q,P)+2j2 A*,(g,P). (16) 

i=2 

© is fully defined once * ~ tpT A, S and are 
specified. S and are primary elements in studies of 
meson properties and are already well constrained. For 
the dressed-quark propagator: 

S(p)=-h- P a v (p 2 )+a s (p 2 ) (17) 
= \i 1 -pA(p 2 ) + B( P 2 )}- 1 , (18) 

we use the algebraic parametrisations |j| : 



2 



a s (x) = 2mT(2{x + rh 2 )) (19) 
+J r (b 1 x) T(b 3 x) [b + b 2 T{ex)} , 
, > 1 



[l -T{2{x + to 2 ))] 



(20) 



with -F(y) = (1 — e~ y )/y, x = p 2 /X 2 , to = to/A, 0-5(2;) = 
X&s(p 2 ) and CTy(x) = A 2 c7y(p 2 ). The mass-scale, A = 
0.566 GeV, and parameter values 



"11 



bi 



0.00897 0.131 2.90 0.603 0.185 



(21) 



were fixed in a least-squares fit to light-meson observ- 
ables. (e = 10~ 4 in (|l^ ) acts only to decouple the large- 
and intermediate-p 2 domains.) This algebraic parametri- 
sation combines the effects of confinement and DCSB 
with free-particle behaviour at large spacelike p 2 B . 

In ([To|)-(|l5|), T M is the dressed-quark-photon vertex. It 
satisfies the vector Ward-Takahashi identity: 



(22) 



r M has been 



which ensures current conservation || 
much studied (l6) and, although its exact form re- 
mains unknown, its qualitative features have been elu- 
cidated so that a phcnomcnologically efficacious Ansatz 
has emerged jL7J: 

iTpVM =m A {£\,tl) lll (23) 
+ (h+l2)» [\i 1 -{h+£ 2 )A A {i 2 ,£ 2 ) + A B (£ 2 1 i 2 )\ ; 

£H«) = i[F(£i) + F(^)], (24) 
F(£ 2 )-F(£ 2 ) 



A F (£(,£ 2 



e 2 



(25) 



where F = A, B; i.e., the scalar functions in (fl7|). A 
feature of ( p3| ) is that r M is completely determined by 
the dressed-quark propagator. Further, we estimate that 
calculable improvements would modify our results by < 
15% Hi. 

The new element herein is the model of the nucleon's 
Fadde'ev amplitude, (||). For the Bethe-Salpeter- like am- 
plitudes we use the one-parameter model forms 



r(gi, ga) = T(q 2 /ul) , Q-= §(<7i -52) 



(26) 



WM = j^F{pl<4), e~{l-ri)h - V £ 2 . (27) 

Our impulse approximation is founded on a dressed- 
ladder kernel in the Fadde'ev equation and satisfies 
(p2|). Hence, the canonical normalisation conditions for 
the diquark and nucleon amplitudes translate to the con- 
straints that the (ud)-diquark must have charge 1 /3 and 
the proton unit charge, which fix Ar and A/*. For the 
diquark propagator we use the one-parameter form 



A(i( 2 ) = i^ 2 /^) 



(28) 



and interpret 1 / toa as the diquark correlation length. 

We fix the model's three parameters by optimising a 
fit to G p E {q 2 ) and ensuring Gjj(0) = 0, which yields^ 



(29) 









TOA 


77 = 2/3 


0.20 


1.0 


0.63 



all in GeV (1/toa = 0.31 fm). Using Monte-Carlo meth- 
ods to evaluate the multi-dimensional integrals, these val- 
ues give 





emp. 


calc. 


r 2 (fm) 2 
rl(im) 2 
M P (Mat) 

Mn//^p 


(0.87)* 
-(0.34) 2 
2.79 
-1.91 
-0.68 


(0.79)* 
-(0.43) 2 
2.88 
-1.58 
-0.55 



(30) 



where the statistical error is < 1 %. The sensitivity of 
our results to the model's parameters is illustrated in 
Table §. It is clear that the fit is stable but does not 
bracket the experimental domain; i.e., the model lacks a 
relevant degree of freedom, a defect we expect including 
an axial-vector diquark to ameliorate. 
The charge radii are obtained via 



p,n 



-6 



Fr(q 2 ) 

\2 



dq 2 

■= (r 1 ) 2 + (r F ) 2 

V p.n/ 1 V p,nJ 

and in this calculation (in fm 2 ) 



2M 2 



if'"(0), 



F\2 



(0.70) 2 , ( r n 

-(0.29) 2 , K) 2 



(0.35) 2 , 
-(0.32) 2 . 



(31) 
(32) 

(33) 



by 



A 20% reduction in loy (Tabic Q, row 4) reduces \r„ 
7%. However, that results from a 21% reduction in |r,_ 
and 2% increase in |r„|. We attribute our overestimate 
of |r 2 | to a poor description of F"(g 2 ), which involves 
many cancellations between terms because of the (u, d, d) 
electric charge combinations and must vanish at q 2 = 0. 

Five diagrams contribute to our impulse approxima- 
tion and diagram 2 involves the diquark form factor. The 
calculated value of the associated elastic charge radius 
provides a measure of the size of the "constituent" di- 
quark: 



r 2 + = (0.45 fm) 2 = (0.80 r^f 



(34) 



3 Our results are sensitive to r\ because ( pq ) and (E7]) are 
equivalent to retaining only the leading Dirac amplitude in the 
expression for these functions and neglecting their q ■ K, I P 
dependence when solving the Bethe-Salpeter and Fadde'ev 
equations. 77 = 2/3 is required for this Ansatz to transform 
correctly under charge conjugation. Accounting for the q ■ K, 
t ■ P dependence would eliminate this artefact [[14|]l9|| . 
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with rv calculated in the same model H, and in quantita- 
tive agreement with another estimate p2j . This is impor- 
tant because, with cur allowed to vary, r + is a qualitative 
prediction of the model. Thus an optimal description of 
the data requires a nonpointlikc diquark. 

Table || provides a guide to each diagram's relative 
importance. In all cases the first diagram, describing 
scattering from the dormant quark, is the most signifi- 
cant. For the charge radii the breakup contributions are 
comparable in magnitude to the second diagram, photon- 
diquark scattering. The magnetic moments are of partic- 
ular interest. A scalar diquark does not have a magnetic 
moment, and that is expressed in our calculation by the 
very small contribution from diagram 2. It is not iden- 
tically zero because of the confinement of the spectator 
quark; i.e., the absence of a mass-shell. Diagrams 3-5 
only appear because the diquark is a nonpointlike com- 
posite and they provide ~ 50% of \i n . Discarding 
these contributions one obtains > —0.5, and in 

pointlike diquark models the axial-vector has alone been 
forced to remedy that defect [^3|. Our results indicate 
that approach to be erroneous, attributing too much im- 
portance to the axial-vector correlation. 

The calculated form factors are depicted in Figs. ^ 
and H and it is obvious in Fig. || that we used G p E (q 2 ) to 
constrain our fit. The + (ud) diquark correlation in 'F 
ensures that G Efit (q 2 ) ^ 0, and the presence of diquark 
correlations can also explain the N-A mass difference. 
Our result for G r ^(q 2 ) is well described by pO] 



G n Efit (q 2 ) 



-u emp F (a 2 

H-n 1 cmp V<i 



a 2 r 
1 + b 2 - 



with r = q 2 /(2M) 2 , F cmp (q 2 ) given in Fig. § and 
a = 1.33, b = 1.00, and the discrepancy between our 
calculation and experiment can be discussed in terms of 
these parameters, a characterises the charge radius and 
it is < 30% too large, as can be anticipated from (30). 
b describes the magnitude at intermediate momenta and 
it is only ~ 23-35% of the empirical value. That is a 
systematic defect shared by other studies Q that only 
retain the scalar diquark correlation. Unlike those stud- 
ies, however, our calculated magnetic form factors, Fig||, 
agree well with the data and, as we have seen, that is be- 
cause we include the diquark breakup diagrams. 

TABLE I. A variation of the model parameters: ui^,, cur 
and 771a (in GeV) illustrates the sensitivity and stability of 
our results. The column labelled "r„" lists: sign(r^) Ir^l 1 ^ 2 . 
(Radii in fm, magnetic moments in units of pjv. The statisti- 
cal errors are < 1%.) 





Tp T n jAp fjj n fj, n f flp 


0.20 1.0 0.63 


0.79 -0.43 2.88 -1.58 -0.55 


0.16 1.0 0.63 
0.24 1.0 0.62 


0.84 -0.46 2.83 -1.55 -0.55 
0.75 -0.41 2.89 -1.59 -0.55 


0.20 0.8 0.62 
0.20 1.2 0.63 


0.80 -0.40 2.93 -1.64 -0.56 
0.78 -0.45 2.84 -1.54 -0.54 
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q . 

FIG. 2. Upper panel: Calculated proton electric form 
factor: +, compared with the empirical dipole fit: 
F cmp {q 2 ) = 1/(1 + 9 2 /™cm P ) 2 , m-cmp = 0.84GeV. Lower 
panel: Calculated neutron electric form factor: +, com- 
pared with the experimental data as extracted using 
the Argonne V18 potential H]. In both calculations the 
Monte-Carlo errors are smaller than the symbols. 

It must be borne in mind that in our calculation a and b 
are not independent. Modifying the parameters in ( pp| ) so 
as to reduce a automatically and substantially increases 
b. However, not withstanding our observation that its 
importance has previously been overestimated, without 
an axial-vector diquark correlation it is not possible to 
accurately describe all observables simultaneously. 

We have employed a three-parameter model of the nu- 
cleon's Fadde'ev amplitude, 'F, to calculate an impulse 
approximation to the electromagnetic form factors. 'F 
represents the nucleon as a bound state of a confined 
quark and confined, nonpointlike scalar diquark, and 
the exchange of roles between the dormant and diquark- 
participant quarks is an integral feature. Five processes 
contribute: direct quark-photon scattering with a spec- 
tator diquark; photon-diquark scattering with a specta- 
tor quark; and three distinct diquark breakup diagrams. 
We obtain a good description of all form factors except 
Gfi, which is too large in magnitude. That defect is 
shared by all models that do not include more than a 
scalar diquark correlation. The nonpointlike nature of 



4 



the diquark correlation is important, especially via the 
breakup contributions which provide large contributions 
to the magnetic moments and ensure fJ, n /^p < —0.5. 
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Energy, Nuclear Physics Division, under contract number 
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A.v. Humboldt foundation. 
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FIG. 3. Calculated proton and neutron magnetic form fac- 
tors, normalised by |/x p ,n| in (|3c|). The curves are dipole fits 
with masses (in GeV): rrip = 0.95, m n — 1.0, 13% and 19% 
larger than m em p in Fig. 0. (/^p™ p F emp (q 2 ) describes the data 
very well.) 

Including a nonpointlike axial-vector diquark is an ob- 
vious improvement of the model. That must be done in 
analogy with the scalar diquark because an accurate in- 
terpretation of the model parameters is impossible if the 
breakup diagrams are discarded. Another avenue for im- 
provement is a direct solution of the Fadde'ev equation, 
retaining the axial-vector correlation and the breakup 
contributions to the form factor. That would provide 
a model for correlating meson and baryon observables in 
terms of very few parameters. 

Models of the nucleon such as ours have hitherto been 
applied only at small- and intermediate-^ 2 . Based on 
the observation || that a description of the large-g 2 be- 
haviour of F^^q 2 ) is only possible if the subleading pseu- 
dovector components of the pion's Bethe-Salpeter am- 
plitude are retained, we anticipate that a successful de- 
scription of the nucleon form factors on that domain will 
require a parametrisation of the Fadde'ev amplitude that 
includes the analogous subleading Dirac components. 

TABLE II. Relative contribution to the charge radii and 
magnetic moments of each of the five diagrams in our impulse 
approximation: Fig. (0), 



diagram 



KY/rl 



1 



0.68 0.11 -0.02 0.12 0.12 

1.14 -0.37 -0.15 0.19 0.19 

0.60 0.01 0.04 0.17 0.17 

0.55 -0.02 0.15 0.16 0.16 
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